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ABSTRACT 

We use our vertical Milky Way disc model together with Galaxia to create mock ob¬ 
servations of stellar samples in the solar neighbourhood. We compare these to the 
corresponding volume complete observational samples of dereddened and binary ac¬ 
counted data from Hipparcos and the Catalogue of Nearby Stars. 

Sampling the likelihood in the parameter space we determine a new hducial initial 
mass function (IMF) considering constraints from dwarf and giant stars. The result¬ 
ing IMF observationahy backed in the range from 0.5 to 10 Mq is a two slope broken 
power law with —1.49 ± 0.08 for the low mass slope, a break at 1.39 ± 0.05 Mq and a 
high mass slope of —3.02 ± 0.06. 

The Besancon group also converging to a similar IMF even though their observational 
sample being quite different to ours shows that the forward modelling technique is a 
powerful diagnostic to test theoretical concepts like the local held star IMF. 

Key words: Galaxy: disc, evolution, formation, solar neighbourhood; Stars: lumi¬ 
nosity function, mass function; Methods: statistical 


1 INTRODUCTION 


The initial mass function (IMF) has seen numerous 


determinations since the seminal work of Salpeter (1955). 


1986 


The most popular being broken power laws (Scalo 
[Kroupa et al.||1993 ) or lognormal with a Salpeter-like high 
mass slope ( Chabrier|2003 ). Originally the IMF was derived 


by luminosity functions (LFs) of the solar neighbourhood 
via the luminosity-mass relation and stellar lifetimes. More 
sophisticated derivations account for the observational sam¬ 
ple’s star formation history (SFH) (|Schmidt|1959|) and con¬ 


nectedly the scale height dilution ( Miller Scalo||l979 ). 

Another more direct approach is to look at young stellar 
clusters and derive an IMF from these spatially and tempo¬ 
rally conhned starburst populations. This yields a range of 
different IMFs ( van den Bergh fc Sher[19^ Dib|2014 ) which 
is not necessarily in conflict ( Weidner Kroupa|2005 ) with 
a global IMF being able to describe the field stars in the 
Milky Way disc. In this paper we try to get hold of the the¬ 
oretical concept of a time independent local Milky Way IMF 
by exploring its effect within our disc model in the realm of 
observables. 


The series of papers started with Just & JahreiB (2010 


hereafter paper I) where a local vertical Milky Way disc 
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model was constructed using dynamical constraints of main 
sequence (MS) stars. The main advantage of this method is 
its weak dependence on the assumed IMF as only the in¬ 
tegrated mass loss goes into the dynamical model. Other 
determinations as for example Haywood et al. (1997) have a 
degeneracy of their derived SFH and their IMF’s high mass 
slope. As a result the old Besangon model (Robin et al. 2003) 
using the constant SFH from Haywood et al. (1997) had 
to be revised recently ( Czekaj et al.||2(3T4 ) to a decreasing 
one in line with different determinations from chemical evo¬ 
lution models (Chiappini et al. 1997), extragalactic trends 


( Ly et aL]|2007 ) and dynamical constraints (|Aumer & Bin- 
ney|2009| [just Jahreifi||2010 ). 


In paper H of this series ( Just et al.|2011 ) the disc model 
was further constrained by comparing Sloan Digital Sky Sur¬ 
vey (SDSS) data with predicted star counts of the north 
galactic pole (NGP). This yielded the fiducial model A (from 
now on called JJ-model) with fixed SFH, a well defined local 
thick disc model (around 6 % of local stars) which showed 
very good agreement between our model and the NGP data 
being confirmed in [Czekaj et al. ( 2014| ). 


After fixing the SFH, age-velocity dispersion relation 
(AVR) and age-metallicity relation (AMR) we now want 
to further build up a consistent disc model by constrain¬ 
ing the IMF using local star counts and taking into account 
turn-off stars and giants as well. In paper I the present day 
stellar mass function (PDMF) was already converted into 
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an IMF considering scale height correction and finite stel¬ 
lar lifetimes. Binarity and reddening were not accounted for 
producing a very steep high mass slope with a = 4.16 which 
is also related to a relatively high break in the power law at 
1.72 M©. 

In this paper we redo the selection of volume complete 


samples from the revised Hipparcos Catalogue (van Leeuwen 


2007) and an updated version of the Catalogue of Nearby 


Stars 4 (JahreiB & Wielen 1997 which will soon be pub¬ 
lished as CNS 5) also giving us a sound statistical sample 
to investigate the local stellar mass density (see discussion 
in sec5.3). We use Galaxia (Sharma et al. 2011) to create 


mock observations from the JJ-model with arbitrary IMFs. 
We find a new fiducial IMF describing the data by sampling 
the two slope broken power law parameter space using a 
Markov chain Monte Carlo technique. The model likelihood 
given the data is approximated assuming discrete Poisson 
probabilities for different magnitude bins differentiating be¬ 
tween dwarfs and giants. 

A similar determination of the IMF was done in the |Czekaj| 


et al. (2014) Besangon update with matching results. In¬ 


stead of local luminosity functions they used the colour pro¬ 
jections of Tycho 2 ( Hpg et al.|2000 ) data in different direc¬ 
tions to determine their best fit. Their modelling machinery 
is very modular being able to incorporate different stellar 
evolutionary- and atmosphere models as well as account¬ 
ing for binarity or reddening from the model’s side. On the 
other hand they have so many degenerate free parameters 
that they do not explore the full IMF parameter space but 
use pre-determined ones from literature. 

Overall the holistic approach with a concordance 
Galaxy model in the background able to produce mock 
observations to test theoretical concepts like for example 
the IMF is a very promising technique especially with the 
underlying physical models getting more and more refined 
with increasing observational evidence. One unsolved prob¬ 
lem which is also pointed out by Czekaj et al. (2014) are 
missing detailed 3d extinction data of the local ISM. 

In sectionj^ we explain how to create mock observa¬ 
tions from our vertical disc model using Galaxia. Sectionj^ 
describes the reduction of Hipparcos and CNS 5 data to ob¬ 
tain our observational sample, followed by sectionj^on the 
likelihood determination in the IMF parameter space. The 
results are then presented in sectionj^ together with a com¬ 
parison to widely used IMFs. In sect ion[^ the results are put 
into context followed by the conclusions in sectionj^ 


2 SYNTHESISING A LOCAL DISC MODEL 

In paper I a self-consistent dynamical vertical Milky Way 
disc model was developed. Combined with constraints from 
star counts of the NGP in paper II a fiducial disc model (JJ- 
model) was chosen. The JJ-model fixes the SFH, the AVR, 
a simple metal enrichment law and from that predicts the 
stellar vertical disc structure in terms of kinematics, star 
counts, ages and metallicities as functions of distance to the 
Galactic plane. 

In order to obtain mock observations of the solar neigh¬ 
bourhood we turn our vertical disc model into a local repre¬ 
sentation used as an input for Galaxia to synthesise stellar 
samples. 
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2.1 The disc model locally 

As we want to compare our JJ-model to volume complete 
stellar samples in the solar neighbourhood we construct 
7 spheres that are disjunct in absolute V magnitude and have 
heliocentric distances from 20 up to 200 pc (see table[^&|^ 
for detailed limits). In the following we explain how we pre¬ 
pare n-body particle representations of our model for each 
sphere such that Galaxia can turn them into mock observa¬ 
tions. 

As our model is evolved over 12 Gyr in 25 Myr steps we 
construct 480 n-body particles (plus one for the thick disc) 
with specific ages, masses and metallicities for each sphere 
separately. These are passed to Galaxia where each particle 
is turned into a simple stellar population (SSP) according to 
a prescribed IMF. The two youngest SSPs are split up age- 
wise so that the youngest stars coming from Galaxia are 
6.25 Myr old. We assign arbitrary phase-space information 
to the particles since the final observable is a local colour 
magnitude diagram (GMD) based on volume complete sub¬ 
samples. 

We will now illustrate the construction of a local n-body 
representation of our model for the 25 pc sphere following 
figure[^ 

The upper panel depicts the global SFH in units of surface 
density per time and the effective scale height, hgff, over age 
for the JJ-model. From that a local mass density, po, for 
each SSP can be calculated which then is multiplied with 
the volume of the 25 pc sphere, V 25 , to obtain 
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age [Gyr] 


Figure 1. The upper panel shows in dashed blue the global SFH 
and in solid red the effective scale height of the JJ-model (at solar 
galactocentric distance). 

In the lower panel the input n-body data which is passed to 
Galaxia to create the 25 pc sphere is displayed. The age, mass 
and metallicity of each SSP are visualised. 


M25pc(G) = V25pcPo(G) 


4 3 SFH(tO- 0.025 Gyr 

-JESS)- 


which is shown in the lower panel of figure[^ as resulting 
n-body particle masses. Colour-coded the chemical enrich¬ 
ment law (AMR) is depicted which is also coming from the 
JJ-model as an analytic function but with an added Gaus¬ 
sian scatter of 0.13 dex standard deviation (tested with GCS 
data, cf. fig. 15 of paper I). 

For the thin disc 3552 M© (A4iMF,thin,25pc) in particle 
mass is passed to Galaxia for the 25 pc sphere 


12 Gyr 


/ ^ou 

M25pc(t)dt — ^ M25pc(tz)- (2) 


OGyr 


This is by construction of the JJ-model the gas mass that 
was used to create the thin disc stars (and in the meanwhile 
also remnants) still residing in the 25 pc sphere. Nowadays 
only a fraction of that is left in stars due to stellar evolution 
(A4pDMF,25pc, cf. section 5.3). In the mass-age distribution 
of the n-body data the peak from the global SFH (dashed 
blue line) around 10 Gyr can still be recognised. The increase 
for younger stellar populations stems from the decreasing ef¬ 
fective scale height confining these stars closer to the Galac¬ 
tic plane (i.e. a bigger fraction of them is found in the local 
sphere, cf. fig. 14 of paper I). 

The thick disc is implemented by inserting a single star- 
burst (i.e. one SSP and equivalently one n-body particle) 
with 6.5 % of the thin disc mass 


AJimF, discs — AllMF,thin + AJimF ,thick — 1.065 • AJimf ,thin (3) 

12 Gyr ago with a metallicity of [Fe/H] = —0.7 resulting in 
a present day thick disc mass fraction of around 5 %. Since 
the local density of the stellar halo is negligible compared 
to the disc we do not consider it as a separate component in 
this work. 

A change in the IMF will affect the mass fraction re¬ 
maining in the stellar component (cf. sec. 2.5 paper I) so we 


introduce a mass factor (mf) scaling our model’s total mass 
(AJimf) to fit the observed stellar mass density 

AJimf = mf • AJimf, discs- (4) 

When going to larger heights above the Galactic plane we 
also need to correct for decreasing vertical density profiles. 
For the 25 pc sphere the deviation from homogeneous den¬ 
sity distribution is still negligible but this changes with 
larger spheres and especially with young stellar populations. 
For example when rescaling the star count in the highest 
magnitude bin in table[^ with a mean age of 0.1 Gyr the 
number must be increased by 54 % to account for the low 
scale height of these stars (cf. figure 2 in paper II). 


2.2 Mock observations with Galaxia 


Galaxia ( Sharma et al.|20'TT ) is a tool to generate mock cata¬ 
logues from either analytic models or n-body data. It already 
has a default Besancon-like model ( Robin et al.||2003 ) but 
the updates ( Robin et al.||2012 Gzekaj et al.||2014 ) are not 
implemented yet. 

In the previous subsection we constructed particles rep¬ 
resenting our model locally which we now pass for each 
sphere separately to Galaxia together with the disjunct 
VMag limits building up a GMD successively. 


So far Galaxia uses Padova isochrones ( Marigo et al. 
|20Q8| ) which have problems reproducing the lower end of the 
MS and the red clump. We include their revised templates 
(PARSEG version 1.2 where only minor differences at 
low mass stars persist (cf. lower MS in figure]^. This re¬ 
maining discrepancy is also pointed out in the release paper 
Ghen et al. ( 2014| fig. A3) but should have negligible effect 


on the star counts in our used magnitude range. 

Binaries, white dwarfs or other remnants are not imple¬ 
mented in Galaxia yet but an update is being planned (pri¬ 
vate communication, Sharma 2015). When inspecting the 
GMD in figurej^a second blue-shifted MS in the synthesised 
catalogue is visible which comes from our distinct thick disc 
metallicity. 

Beside being able to change the IMF according to which 
Galaxia is distributing the particle masses into stars we are 
using it as a black box. Specifying a photometric system 
will already yield a detailed stellar catalogu^in terms of a 
random realisation of our local model representation as for 
example depicted in figurej^for our newly determined IMF. 


3 OBSERVATIONS 

The anchoring point for every Galaxy model in terms of ob¬ 
servational constraints is the stellar distribution in the solar 
neighbourhood since detailed and volume complete samples 
can only be obtained here. After using the vertical com¬ 
ponent of the velocity distribution of MS stars in paper I 
(dynamical constraint) and the NGP star counts in paper II 
(vertical density distribution constraint) we now use abso¬ 
lute local stellar densities for dwarf and giant stars. 


^ http://stev.oapd.inaf.it/cgi-bin/cmd 

^ See http://galaxia.sourceforge.net/Galaxia3pub.html for de¬ 
tailed instructions 
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Table 1. Observational sample and mock catalogues - dwarf stars 


Catalogue 

d 

(pc) 

Mv^-limits 

(Mag) 

# 

^ > 15% 

# lost 

CNS 5 
25 pc 

N 25 JJ25 

rescaledi to 25 pc 

log-likelihood 
In (£/Pmax) 

Mean mass 
Mo 

Mean age 
Cyr 


200 

1—1 

1 

98 

0 

0 

0.28 

0.13 

-21.9 

6.4 

0.1 


200 

[-1.5,-0.5] 

233 

3 

1 

0.62 

0.63 

0.00 

4.0 

0.2 


200 

[-0.5,0.5] 

901 

12 

4 

2.26 

2.60 

-9.20 

2.9 

0.3 

Hipparcos 

100 

[0.5,1.5] 

520 

2 

15 

8.61 

9.09 

-0.76 

2.2 

0.5 


75 

[1.5, 2.5] 

677 

1 

27 

25.6 

24.4 

-0.81 

1.7 

1.0 


50 

[2.5, 3.5] 

518 

1 

62 

65.1 

59.9 

-1.79 

1.4 

2.3 


30 

[3.5,4.5] 

200 

0 

110 

115.7 

146.5 

-5.90 

1.1 

5.3 


25 

[4.5, 5.5] 

191 

4 

191 

191 

190.3 

0.00 

0.9 

6.1 


25 

[5.5, 6.5] 

198 

11 

198 

198 

207.9 

-0.22 

0.8 

6.4 

CNS 5 

25 

[6.5, 7.5] 

193 

16 

193 

193 

196.0 

-0.02 

0.7 

6.5 


25 

[7.5, 8.5] 

207 

13 

207 

207 

205.3 

-0.01 

0.6 

6.6 


20 

[8.5, 9.5] 

139 

15 

245* 

271.5 

258.0 

-0.20 

0.5 

6.7 


E 

],9.5] 

4075 

78 

1253 

1278.8 

1300.8 

-40.8 

0.8 

6.0 


Catalogue shows from which source the observational sample is drawn, d gives the heliocentric distance of stars included (sphere 
size), My-limits gives the magnitude range of each bin, Nfin is the final star count in each bin, ^ > 15 % is the number of stars 
thrown out due to high distance errors, CNS5 gives the number of stars within the volume complete 25 pc sphere, the next two 
columns give the star counts of the observations and our JJ-model rescaled to 25 pc. The JJ-model with newly determined IMF is 
averaged over 400 random realisations, log-likelihood shows the probability of each bin after equation]^ normed with the maximal 
possible probability (cf. section [TT] > in natural logarithm which indicates each bin’s impact on the likelihood function (we will call it 
penalty because our likelihood can be seen as a reward function for the MCMC simulation), the mean mass and mean age show the 
values for the corresponding JJ magnitude bins where the sum at the bottom is an average of all stars within 25 pc. 

*for this magnitude bin volume completeness is not given so the 139 stars from the 20 pc sphere have been rescaled to 25 pc yielding 
271.5 stars which is 10% more than the 245 stars observed in the 25 pc sphere 

%escaling the volume and accounting for the density profile of the magnitude bin’s mean age population 


Table 2. Observational sample and mock catalogues - giant stars 


Catalogue 

d 

(pc) 

Mv^-limits 

(Mag) 

^ fin 
# 

^ > 15% 

# lost 

CNS 5 
25 pc 

N 25 JJ 25 

rescaledi to 25 pc 

log-likelihood 
In (C/Pmax) 

Mean mass 
M© 

Mean age 
Cyr 


200 

1,-1.5] 

74 

1 

0 

0.16 

0.14 

-0.87 

3.4 

1.6 


200 

[-1.5,-0.5] 

375 

2 

1 

0.77 

0.84 

-1.47 

1.8 

4.4 


200 

[-0.5,0.5] 

1341 

20 

3 

2.78 

2.29 

-23.4 

1.7 

3.9 

Hipparcos 

100 

[0.5,1.5] 

526 

0 

6 

8.33 

9.27 

-3.03 

1.5 

4.8 


75 

[1.5, 2.5] 

126 

0 

5 

4.70 

3.76 

-2.99 

1.2 

6.3 


50 

[2.5, 3.5] 

62 

0 

7 

7.77 

11.6 

-5.55 

1.2 

6.6 


30 

[3.5,4.5] 

9 

0 

6 

5.21 

7.11 

-0.34 

1.0 

8.8 


E 

],9.5] 

2513 

23 

28 

29.7 

35.0 

-37.7 

1.3 

6.3 


This is achieved by constructing different samples combining 
absolute magnitude cuts and heliocentric distances such that 
the selected stars represent a volume complete sphere. At the 
bright end we go up to 200 pc in order to obtain enough mas¬ 
sive stars and giants to have reliable statistics. Our observa¬ 
tional sample consists of stars from the extended Hipparcos 
catalogue ( [Anderson fc Francis||2012 ) and the Catalogue of 
Nearby Stars (CNS). CNS 5 which we use here is an updated 
version of CNS 4 ( jJahreifi Wielen|1997 ) and will get pub¬ 
lished soon. 

Fundamentally we would like to implement all observational 
biases on the models side and compare the synthesised mock 
observations to unaltered observables. In this respect the up¬ 
dated Besangon model (Czekaj et al. has pushed the 

link between model and data in the right direction by imple¬ 
menting extinction models and a scheme for binary systems 
into their model. As Galaxia is not able to account for bina¬ 
ries yet and we are not providing positional information with 


the n-body particles we have to treat binaries and deredden¬ 
ing from the observational side. 


3.1 Hipparcos 

We use the extended Hipparcos compilation (117955 entries) 
which cross-matches the original stars from the revised Hip¬ 
parcos catalogue ( van Leeuwen||2007 ) with a large selection 
of different catalogues. 

By using heliocentric distance and absolute magnitude cuts 
similar to paper I ( Just Jahreifi|[2010 table 1) we obtain 
volume complete observational spheres for different stellar 
magnitudes going down to 4.5VMag. The only further se¬ 
lection criteria to eliminate miss-identifications is a distance 
error below 15% which reduces the sample insignificantly as 
visible in table[T]&[2l 

Before the VMag and distance cuts are applied all stars 
which have both a binary flag and a Amag entry in the 
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Table 3. Effect of binarity and dereddening on the star counts 


Radius of sphere [pc] 
Magnitude bin [VMag] 

30 

4th 

50 

3rd 

75 

2 nd 

100 

1 st 

200 

< 0.5 

No correction 

204 

591 

793 

1052 

2798 

Binary correction 

209 

580 

803 

1046 

2756 

Extinction correction 

204 

591 

793 

1053 

3060 

Both corrections (Nj^^) 

209 

580 

803 

1046 

3022 


original Hipparcos catalogue (ESA 1997) are split up into 
two components. The binary correction changes the number 
of stars in each distance bin slightly as can be seen in tablej^ 
With the binary correction we have 5394 stars in the Hip¬ 
parcos sample of which 552 entries come from split up binary 
systems. On average the bins lose stars when correcting for 
binaries because the single components get fainter than the 
system bringing these stars below the magnitude limit. As 
the fainter mag bins have smaller limiting radii only some 
lost stars fall into the fainter magnitude bin. 

On scales of the investigated volume the ISM is in- 
homogeneously distributed with Ophiuchus and Taurus 
molecular cloud being the biggest absorbers in the 200 pc 
sphere ( Schlahy et al.||2014 ). 3d extinction maps with good 


enough resolution are getting published (Lallement et al. 


2014) but the data are not available yet. In order to dered- 


den our stars we adopt an analytic model from |Vergely et ah] 
(1998) describing a homogeneous extinction depending on 
the distance and Galactic latitude. Due to the local bubble 
we set the extinction to 0 below 70 pc distance and above 
52° Galactic latitude (c f. [Vergely et al. (1998 hg. 4, 11)). 


Aumer & Binney (2009) also adopted this model but 


only within 40 pc of the Galactic mid-plane resulting in 4 % 
less stars in the 200 pc sphere compared to our adaptation 


of the extinction model. We use the Vergely et al. (1998 
p.548) cosecant law: 


{ 0, d < do V 5 > 52° 

Eo(d-do), if (5) 

Eo (|sin( 6 )| - do) , d> | 3 in( 6 )| 

with ho = 55 pc, do = 70 pc, Eo = 4.7 x 10“^mag/pc 
and d, h being the heliocentric distance and the Galactic 
latitude. To transform from E {B — V) to A (V) we adopt 
Ro = 3.1. 

In essence the dereddening leaves the closer samples unal¬ 
tered only increasing the 200 pc sample star count by around 
10 %. This is due to stars with heliocentric distance between 
100 and 200 pc that are slightly too faint without deredden¬ 
ing become a bit brighter and thereby satisfy the magnitude 
limits. 

As the ISM is highly inhomogeneous adopting an analytic 
model is only a crude approximation but the best we can do 
at the moment and we are curious to redo the analysis with 
upcoming 3d extinction maps. 
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Figure 2. CMD of the observational sample (Hipparcos stel¬ 
lar systems are not split up) and one random realisation from 
our model with the newly determined IMF. The cuts with con¬ 
nected triangles at (B-V,My) = (-0.3,-5.5), (0.5,0.5), (0.75,4.5) 
and (3,4.5) for the division into dwarf and giant sample are in¬ 
dicated in magenta. In the right the different sphere radii are 
written in units of pc. 


The distance of the stars is calculated from parallaxes (when 
advisable the photometric were preferred to the trigonomet¬ 
ric) with a correction for the Lutz-Kelker bias from Ander-| 


[son V Erancis (2012 eq. 1). 

The 20 pc sample includes stars from 8.5 to 9.5 VMag. After 
excluding 2WDs and 15 stars with relative distance errors 
above 15% we are left with 139 stars. Eor the 25 pc sphere 
the absolute V magnitude ranges from 4.5 to 8.5. Here 44 
stars have a too high error so that 789 stars remain in the 
hnal sample. 

Ah together 928 stars come from the GNS5 of which 372 
are in resolved multiple stellar systems (24 of the remaining 
556 single systems are detected spectroscopic binaries). Of 
the 372 stars in multiple stellar systems 117 mostly primary 
and 87 mostly secondary components share a joint B-V value 
but have individual VMag entries so that they can be cor¬ 
rected by putting them on the MS. The MS was empirically 
assigned from MS stars with low parallax error. 

In the end three stars reside in the GNS5 sample as 
well as in the Hipparcos sample because of slightly different 
VMag values directly at the VMag borders of 4.5 VMag. 
Those stars are excluded from the Hipparcos sample so that 
in our joint sample every star has a unique entry. 

As a comparison to the Hipparcos sample we also list the 
247 volume complete (within 25 pc) GNS5 stars which are 
brighter than 4.5 VMag in table[^and|^ The only peculiarity 
is the 2cr outlier with 15 against 8.65 expected stars in the 
100 pc bin of the dwarf sample. Still the GNS 5 sample can be 
seen as a valid random realisation of the enlarged Hipparcos 
sample. 


3.2 Catalogue of Nearby Stars 5 

The Hipparcos sample is supplemented with stars from 
GNS 5 (7251 entries) at the faint end (4.5 - 9.5 VMag) which 
is a volume complete catalogue for stars brighter than 8.5 
(9.5) VMag up to a distance of 25 (20) pc. 


4 STATISTICAL ANALYSIS 

Eor our analysis we divide the derived GMDs into dwarf and 
giant stars with the cuts specihed in hgurej^ Since MS stars 
have a tight correlation between luminosity and stellar mass 
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the dwarf sample contains information on the PDMF. The 
giant sample adds constraints for the integrated SFH and 
the IMF of higher mass stars (above 0.9 M©). 

To compare our observable (the LF) to a theoretical IMF 
we feed local representations of our model together with 
different IMFs to Galaxia leaving us with synthesised star 
counts from which we construct a likelihood assuming dis¬ 
crete Poisson processes (section [4T] ) . This is implemented 
into a Markov chain Monte Carlo (MCMC) scheme to ob¬ 
tain a representation of the probability distribution function 
(PDF) in the two slope IMF parameter space (section [T^ . 


4.1 Likelihood calculation 


We approximate the likelihood of our model given the data 
by dividing each CMD into 12 magnitude bins for the dwarf 
sample and 7 for the giants (see table[^&|^ and calculate 
the discrete Poisson probability distribution. The expected 
value is coming from our model (rrii) and the number of oc¬ 
currences is the star count observed in each bin (di) leading 
to the likelihood 

12+7 di -mi 

Aotai = ri where A = ^ -. (6) 

1=1 

The log-likelihood then follows as: 


12+7 

log£totai = ^ (di \og(mi) - log(dd) - (7) 

i = l 


For the calculation of the log factorial a very accurate ap¬ 
proximation for n > 0 from Aiyangar ( |1988 p. 339) is used: 


logn! nlogn —n + 


log + 4n(l + 2n)^^ 


log(7 


We will normalise the log-likelihood with its maximal possi¬ 
ble value, Pmax = —68.5, occurring when the observed sam¬ 
ple is tested with itself. 

It should be kept in mind that a linear increase in star counts 
results in exponentially increasing penalties for our likeli¬ 
hood function when the deviation in percent stays the same. 
For example if the expected value is 10 and the number of 
occurrences is 9 then we have a ln(£/Pmax) of —0.04. For 
100 expected stars and 90 occurrences ln(£/Pmax) equals 
—0.46, for 1000 and 900 it is —5.1 and so forth. This on the 
one hand takes into account that bins with a lot of stars 
get a higher statistical weight but could also be dangerous 
when small errors in rescaling (which could come from a bad 
AVR or connectedly isochrones indicating wrong ages) result 
in large penalties for the likelihood potentially pointing our 
MCMC simulation to a biased equilibrium IMF parameter 
configuration. 


reduced by a factor of 20 (cf. tab|^ and should be a second 
order effect compared to the Poisson noise in the data (if 
we assume the observed stars have been randomly realised 
from an underlying probability distribution). 


To sample the PDF of our parameter space we use a 
Python implementation (Foreman-Mackey et al. |2013 ) of 


an affine invariant ensemble sampler for MCMC (Goodman 
Weare]|2010 ) where step proposals using the information 


of multiple walkers reduce the autocorrelation time signifi¬ 
cantly. 

Since the overall mass turned into stars (A4imf, discs) could 
be slightly different to the JJ-model we add the mass factor 
(mf) as a fourth free parameter beside the three two slope 
IMF parameters, low mass index (ai), high mass index 
((+ 2 ), and the power law break (mi): 


dn / <+ = (+1, .. miow < m < mi . . 

^ = kam < if (9) 

dm [ a = a2, mi < m < mup 

with the lower and upper mass limit of the IMF miow = 
0.08 M© and mup = 100 M© being fixed and ka normalising 
the IMF to be a continuous function that represents the 
mass turned into stars 

AtiMF = mf • AdiMF, discs = / kam dm. (10) 

-'miow 

For every set of parameters we use the product of the like¬ 
lihoods from dwarf and giant sample 

log/Itotal — lo§-^dwarf T log/^giant (H) 


to sample the parameter space. In this way the probability 
for each bin (12 from the dwarf sample and 7 from the gi¬ 
ant sample) being an independent discrete Poisson process 
is weighted equally into the final likelihood (£totai)- 
We also investigated three slope IMFs but the gain in like¬ 
lihood is small. Usually a second break in the IMF is in¬ 


troduced to fit the low mass regime (Kroupa et al. 1993 


m < 0.5 M©). With our observational evidence ranging from 
0.5 to 10 M© a second break in the IMF below this limit can 
not be tested. But an empirical extension of our 2 slope IMF 
to fulfil local mass density constraints for low mass stars will 
be introduced in sectionlffH 


5 RESULTS 

Here we present the newly determined fiducial IMF for our 
model. For comparison we also show the log-likelihood of the 
observational data with synthesised data generated from our 
model but using common IMFs from the literature. 


4.2 Sampling the likelihood distribution 

The variability of the outcome of Galaxia is twofold. First 
the IMF parameters can be varied changing the laws 
according to which the mock observations are produced. 
Second for fixed parameters the random seed of Galaxia 
can be changed yielding different random realisations. The 
latter can be minimised by averaging. We oversample each 
point in parameter space 400 times so that this noise is 


5.1 New IMF parameters 

Multiple burn-ins from different starting points all settling 
in the same equilibrium configuration suggest a well-behaved 
parameter space with respect to our likelihood calcula¬ 
tion. We use 20 walkers each sampling 500 steps to gener¬ 
ate the point cloud representing the PDF of the parame¬ 
ter space. For each step we average over 400 realisations. 
Figurej^ shows the marginalised likelihood distribution of 
each parameter as histograms on the diagonal and as point 
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a.2 nil mf 


Figure 3. Marginalised parameter distribution of the MCMC 
run exploring the equilibrium distribution of the parameter space 
with respect to our log-likelihood using 10^ evaluations. Each 
scatter plot shows the projected 2d parameter distribution with 
points coloured by likelihood increasing from blue to red. Crosses 
indicate the mean values and ellipses encompass the 1-3 cr regions. 
The respective correlation coefficients are given at the position 
mirrored along the diagonal. Gaussian fits and histograms of the 
marginalised parameter distribution are given on the diagonal. 
The mean and standard deviation of each parameter is written 
and also indicated by solid and dashed grey lines. 





Figure 4. Luminosity function of the observations and the new 
IMF of the JJ model (cf. table[^and[^. Error bars indicate Pois¬ 
son noise in the observational sample and the standard deviation 
for the 400 times oversampled synthesised catalogue. Star counts 
are not normalised for the different distance limits. The limiting 
radii of the corresponding magnitude bins are written in the top. 


clouds for each parameter pair in the lower left. Just for 
illustrative reasons each dot is coloured to indicate high 
(redmax = —77.2) and low (bluemin = —86.7) log-likelihood 
values with grey dots being below this 3 a range. In the his¬ 
tograms the mean and the standard deviation of the number 
density of each marginalised parameter is given which rep¬ 
resents our central result and dehnes our new hducial IMF: 

ai = 1.49 zb 0.08, 

0^2 = 3.02 zb 0.06, . . 

mi = 1.39 zb 0.05, ^ ^ 

mf = 1.09 zb 0.04. 


Furthermore Pearson’s correlation coefficient for each pa¬ 
rameter pair is given in the upper right of hgurej^ which is 
also represented in the 1-3 cr ellipses of the projected point 
clouds. The only two parameter which are uncorrelated are 
the power law indices. Almost positive linear is the corre¬ 
lation of the mass factor with the low mass slope. This is 
due to more mass being put into stars which are not rep¬ 
resented in our data (m < 0.5M©) for high ai which can 
be counterbalanced with a high mass factor. Similarly but 
less strong the anti-correlation of the mass factor with the 
high mass slope is due to mass being shifted out of our rep¬ 
resented data domain when a 2 is getting lower. The reason 
for both power law indices being positively correlated with 
the power law break is due to the shape of the IMF which 
is sharply decreasing at the position of the break and the 
number densities of the faintest and brightest magnitude bin 
which need to be matched when shifting mi. This reasoning 
can be best visualised when looking at the green dotted line 
(our IMF) and the blue error bars (our observational sam¬ 
ple) in hgure[^ Additionally when increasing both ai and 
0^2 the mass that is gained by the steeper low mass slope 
will be counterbalanced by the mass loss of the steeper high 
mass slope. 

Figurej^ shows the binned luminosity function of our hdu¬ 
cial IMF (averaged over 400 realisations) compared to the 
observations in absolute numbers. The CMD representation 
of the luminosity functions can be seen in hgurej^ where the 
same colours have been used for the different samples. The 
deviations in each bin look small and systematics are not 
apparent neither in the dwarf nor the giant sample which 
shows that the whole machinery, consisting of the disc model 
and Galaxia producing the mock observations, works well 
and the MCMC simulation has likely converged towards the 
equilibrium conhguration. 

When inspecting the detailed likelihood contribution of each 
bin in table[^&|^we see that the largest penalty comes from 
the 0th VMag bin. Especially the giants with In (£/Pmax) = 
—23.4 have a huge impact. The reason for that is likely a red 
clump that is too faint in our mock catalogue. In the dwarf 
sample too many stars are in the synthesised 0th VMag bin 
with JJ 25 -(N/in/N 25 ) = 1036 compared to 901 stars in the 
Hipparcos catalogue. 

A weakness of our modelling machinery is apparent in the 
high In (£/Pmax) value of the brightest dwarf bin indicating 
that too few bright stars are produced. The reason is most 
likely that we are not accounting for binaries. Minor effects 
could be missing high metallicity stars and that our syn¬ 
thesised stars are not younger than 6.25 Myr. It could also 
indicate a change of the high mass power law index for stars 
more massive than contained in our data. 

5.2 Tested IMFs 

Because for other shapes of the IMF part of the mass could 
also be hidden in the mass range not represented by our 
observational sample (0.5 to 10 Mq) we determine the fac¬ 
tor with which the total mass needs to be rescaled in order 
to maximise the likelihood for a particular IMF when it is 
used with our JJ-model against the observational sample. 
We again average over 400 realisations and get the standard 
deviation as the enclosing 68 % of the likelihood. 

In table[^the log-likelihoods of the different IMFs (using our 
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Table 4. Likelihoods of the different IMFs 

IMF JJ JJsa BesB KTG 93 Cha03 

-79.3 -86.7 -96.7 -195.6 -216.1 

\ ^max / 


Table 5. Variability of the log-likelihood 

Sample JJ 400 JJl JJideal* 

In f -79.3 ±0.6 -89.5 ± 14 -13.8 ±2.8 

*if data was coming from our model (see text) 


disc model but adjusting for each IMF’s best fit mass factor) 
are listed with JJsa being the log-likelihood value which is 
lower than 99.7% of the points representing the PDF in hg- 
ure[3l 

Table[^ illustrates a few properties of our log-likelihood. 
JJ400 and JJl show the mean and the standard deviation of 
100 log-likelihood determinations with different seeds which 
are averaged over 400 in the former and 1 realisation in the 
latter case. This shows that the averaging is important for 
the MCMC simulation in order to smoothen the likelihood 
distribution. The deteriorated mean for single realisations is 
due to a skewed distribution since Pmax is a lower limit and 
the increasing penalty for extreme values. 

The last column of table[^ (JJideai) gives an ideal log- 
likelihood which is obtained when the data is indeed rep¬ 
resented by the model. For that we draw 100 single ran¬ 
dom samples (JJi) from JJ400 with replacement and evalu¬ 
ate their log-likelihood with the parent distribution (JJ 4 oo)- 
Each sample fulfilling the observational constraint of having 
dwarf and giant star counts hxed to Nobs,dwarf = 4075 and 
Nobs,giant = 2513. This means that a In (£/Pmax) of around 
— 13.8 would indicate a perfect model. An even lower log- 
likelihood value close to Pmax (ln(£/Pmax) = 0) would be 
unrealistic since there is a natural scatter to Poisson pro¬ 
cesses. 

Related to that we also inspected the distribution of star 
counts in individual magnitude bins for random realisations 
(with the same parameters but different seeds) which indeed 
is Poissonian. 


KTG 93 


The widely used Kroupa et al. (1993) IMF is a three slope 
broken power law of the form 


( a = ai, miow <'^^ < mi 

-— = kam~^ < a = a2, if mi < m < m 2 (13) 

dm 

a = 0^3, m 2 < m < mup 

with ai = 1.3, a 2 = 2.2, as = 2.7, mi = 0.5 and m 2 = 1. 
Again is a normalisation constant ensuring a continuous 
IMF between miow = 0.08 M© and mup = 100 M©. 

The likelihood peak is obtained with a mass factor of 1.392± 
0.019 which is quite high. The reason for that is too much 
mass being put into low mass stars which are not represented 
in our observational sample (see table|^. 

With ln(£/Pmax) = —195.6 it scores poorly compared to 
our or the BesanconB IMF showing that its shape is not 
able to reproduce local star counts. This is visible in hgure[^ 


where it produces too few stars within 0.9 - 2M© and too 
many outside of this range compared to our IMF. 


BesangonB 

This is one of the new hducial Besancon model IMFs from 


Czekaj et al. (2014) tested with Tycho 2 all-sky colour dis¬ 


tribution. It is also a three slope broken power law with 
ai — 1.3, a 2 = 1.8, as = 3.2, mi = 0.5 and m 2 = 1.53. 

The mass factor for this IMF is 1.107±0.015 which is com¬ 
patible with our own mf. Also the shape of the IMF (cf. 
hgure[^, the likelihood (see table[^ and the mass fractions 
(cf. tablej^ are similar. This is not surprising since they 
are also using the local mass density model from [Jahreifi fc] 
Wielen| ( [l997| . Still it is reassuring that for different observa¬ 
tional data (Tycho 2 colour vs. Hipparcos/CNS5 VMag) and 
different modelling techniques similar results are yielded. Al¬ 
though the difference of 17.4 in In (£/Pmax) still means that 
the likelihood for their IMF is 36 million times lower than 
for our IMF. 

Compared to our IMF the BesanconB IMF is producing 
slightly less high mass stars and more low mass stars which 
could be partly due to their rigorous treatment of binaries 
(see hgurej^and cf. section 6.3). 


Chahrier 03 


Another widely used IMF comes from Chabrier (2003). It is 
a mixture of a lognormal form in the low mass and a Salpeter 
power law in the high mass regime 


dn 

dm 


exp 


0.079) 

2-0.692 


0.237912 • m“^ % 


if 


m < M© 
m > M©. 


(14) 


The mass factor for the best likelihood is: 1.317 ±0.017. It 
is so high because a huge mass fraction is going into stars 
more massive than 8M© (cf. table[^ highest supernova 
rate compared to any other standard IMF). The shape 
of the Chabrier 03 IMF hts our model worst with respect 
to the data scoring a log-likelihood of In (£/Pmax) = —216.1. 


5.3 Prom luminosity function to local stellar mass 
density 

In hgurej^ the IMFs are normalised such that their inte¬ 
grated mass is representing the mass of gas that was turned 
into stars (of which a few already turned into remnants) 
still residing in the 25 pc sphere (i.e. A4iMF,25pc including 
the thick disc and each IMF’s mass factor cf. equation]^. Be¬ 
ware that in hgurej^ actually the number of stars per mass 
interval is displayed though normalisation happens in mass 
space. 

The yellow shaded area shows the incompleteness of our ob¬ 
servations in low mass stars due to the magnitude cut. This 
comes from high metallicity stars being redder and fainter 
than their equal mass low metallicity counterparts so they 
are excluded and do not contribute to the smallest mass bin 
leading to a decrease of the observed mass compared to the 
theoretical IMF. The same is true for the JJ mock dwarfs 
in dotted red (as we apply the same magnitude cuts) which 
were spawned using the green dashed IMF. We could drop 
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Figure 5. 25 pc LF of dwarf stars translated into mass space. 
Our fiducial IMF is plotted in thick dashed green and the Be- 
sangonB and KTG93 IMFs are plotted in dotted and dashed 
grey. The blue error bars represent the 12 magnitude bins from 
the observational dwarf sample with their limiting radius given in 
parsec. The x-value and the x-error are associated with the me¬ 
dian and the range of stellar masses in the corresponding mock 
magnitude bin synthesised with our fiducial IMF. The y-value 
represents the number of stars normalised to the mass range and 
the y-error is the Poisson noise. As the masses from different mag¬ 
nitude bins overlap the values are added up in the solid cyan line 
representing a kind of observational PDMF. Below 0.7 Mq the 
yellow shaded area indicates the mass that is missing due to the 
magnitude cut. In dotted red the same effect is visible as this 
line represents the synthesised mass function of our IMF aver¬ 
aged over 100 realisations. From 0.9 Mq upwards the IMF and 
the local IMF in solid blue deviate since locally (i.e. close to the 
Galactic plane) young and therefore massive stars are over repre¬ 
sented which is indicated by the blue shaded area. Therefore the 
PDMF in the 25 pc sphere would look like the local IMF if there 
was no stellar evolution (shown in magenta) at work. 


the magnitude cut for faiut stars iu Galaxia aud the red dot¬ 
ted hue would be perfectly aligued with the JJ IMF iu the 
low mass regime which correspouds to the local preseut day 
mass fuuctiou (PDMF). Iu the observational PDMF repre- 
seuted by the solid cyau hue the little bump compared to 
the JJ mock dwarfs at arouud 0.8 Mq could be related to 
the isochroues badly fittiug our low mass stars biasiug the 
associated mass rauges. 

Iu the high mass regime the blue shaded area shows the 
over represeutatiou of massive stars iu the local IMF of the 
25 pc sphere as their vertical distributiou is more coufiued 
to the Galactic plaue (if oue would uot accouut for the scale 
height dilutiou the deduced IMF would look like this). The 
mageuta area iudicates the deviatiou from the local IMF to 
the mock sample due to stellar evolutiou. 

We cau quautify the differeuce of IMF aud PDMF by look- 
iug at the iutegrated mass of the two fuuctious represeuted 
by JJ IMF aud JJ mock dwarfs from figure]^ (the latter 
beiug made equal to JJ IMF for the low mass part uot 
affected by stellar evolutiou). The outcome is that 55.8% 
of the mass origiually turued iuto stars is still preseut iu 
dwarf stars today (A4pdmf,25pc/AJimf,25pc )• lucludiug gi- 
auts this value ouly chauges slightly to 56.9 % which is cou- 
sisteut with the iutegrated stellar mass of our disc model 


(figure 1 of paper I adapted for 25 pc age distributiou yields 
56.5% of stars iucreasiug to 70.8% wheu remuauts are iu- 
cluded which is their geff) coufirmiug that the uew IMF still 
fits withiu our model’s framework as SFH aud AVR are ouly 
depeudeut ou the iutegrated mass loss. 

Iu order to clarify our remuauts which Galaxia is uot 
syuthesisiug we derive from our PDMF aud IMF that we 
should have 707 stellar remuauts for stars betweeu 1 aud 
8Mq (poteutially WDs) aud 17 heavier oues (poteutially 
black holes or ueutrou stars) iu the 25 pc sphere. Gompared 
to Siou et al. (2014) we have over a factor of 2 more siuce 


they expect 344 WDs withiu the same limits. But their as¬ 
sumed volume completeuess for the 13 pc sample seems quite 
optimistic at least wheu speakiug of the cool eud of the WD 
cooliug sequeuce. |Holberg et al. (2008) propose a WD meau 
mass of 0.665 Mq which for us result iu a WD mass deusity 
of 0.0072 Mq/pc^ quite close to the Besaugou value based 
ou |Wieleu| (|1974|) wh ich was already corrected dowuwards 
iu [Jahreifi Wieleu (1997) because oue out of 5WD left 
the 5 pc sphere. We admit that our high uumber of WDs is 
partly due to the two slope IMF structure haviug a bump 
where the power law break lies (wheu expectiug a coucave 
fuuctiou describiug the uuderlyiug distributiou iu log-space) 
which slightly exaggerates the mass fractiou of the IMF go- 
iug iuto plauetary uebulae (cf. table[^. 

With these values (summarised iu table[^ the over¬ 
all preseut day mass fractiou of stars aud stellar remuauts 

(WlpDMF,total, 25pc) IS 


5'eff = 


AlpDMF, 25 pc + AJbD &WD 25 pc 


M 


= 71.7% 


IMF,25 pc 


(15) 


which is close to 70.8% for the 25 pc sample of paper I. 
This also implies a combiued BD aud WD mass fractiou 
of about 20% of the local mass budget (A4pDMF,totai, 25pc) 
cousisteut with the Jahreifi & Wieleu (1997 tab. 3) value 
aud the stellar evolutiou of paper I. Theu agaiu we could 
als o have choseu the p roposed WD local mass deusity 


of Holberg et al. (2008) Mq/pc*^ = 0.0032 which would 
decrease the local disc mass deusity aud chauge geff aud the 
remuaut fractiou. 

For the followiug local mass deusity test we use our disc 
model aud the fiducial IMF to produce all stars withiu 
25 pc dowu to 0.08 Mq without auy maguitude cuts aud 
aualyse the sample’s properties. Iu this sample we fiud 
a local stellar mass deusity of 0.034 Mq/pc^ co mpared 
to 0.030 Mq/pc^ for the same selectiou of stars iu Flyuu 
et al. ( |2006| tab. 2) (iucludiug giauts, excludiug BDs, WDs, 
other remuauts aud the thick disc compoueut) who use a 
similar method. Iu table[^ the detailed comparisou reveals 
that especially the two brightest bius of the FH06 sample 
are 3 times deuser thau the correspoudiug bius iu the JJ 
sample. For VMag > 3 the deusities are matched quite well 
except for the faiutest mag biu. The over abuudauce of 
bright stars could be comiug from their large value of u iu 


Holmberg et al. (1997 eq. 3) resultiug iu uearly expoueutial 
vertical deusity distributiou with high local mass deusities 
though this theu should apply to the other mag bius as 
well. Auother iudicatiou for au over estimatiou of their 
model’s local star couuts becomes evideut wheu compariug 
the uumber deusities from Holmberg & Flyuu ( 2000| (upou 
which FH06 is based) with the GNS5 oues. For My < 2.5 
they have 0.0013 star pc“^ whereas our volume complete 
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Table 6. Local stellar mass density of different thin disc compo¬ 
nents in 10 “‘^Mq/pc^ 


Mass component 

Flynn et al. 

(2006|) JJ-model 

giants 

6 

7 

Mv < 2.5 

31 

11 

2.5 < Mv < 3 

15 

5 

3 < Mv < 4 

20 

19 

4 < Mv < 5 

22 

26 

5 < Mv < 8 

70 

66 

8 < Mv 

135 

205 

Pthin disc (-2^ = 0 pC, t = 12 Gyr) 

299 

338 

Pthickdisc(2; = 0pc,t = 12 Gyr) 

35 

17 

brown dwarfs 

20 

20 * 

white dwarfs 

60 

72t 


Ppdmf, total 414 447 

* taking the same value as FH 06 
^derived implicitly from our PDMF (see text) 


sample has half of it with 0.0007 star pc“^ and this despite 
the fact that the 25 pc from the CNS 5 seem to have an 
over representation of the upper MS compared to the larger 
sample (cf. table[^. The next mag bin 2.5 < My < 3 is 
three times denser according to Holmberg & Flynn (2000) 
with 0.0010 star pc“^ compared to 0.0003 star pc“^ we 
measure for our 25 pc sample which strongly indicates a 
necessary revision. Fainter mag bins are much better fit in 
star counts as well as in stellar mass density. 

Comparing our thin disc stellar mass density (excluding 
thick disc stars and WDs) of 0.0338 Mq/pc^ to the default 
BesangonB model that has 0.0330 Mq/pc^ reveals a sim¬ 
ilar discrepancy. Our disc model has 12 Gyr of evolution 
compared to 10 Gyr and we make a detailed comparison 
in table[^ keeping in mind the different SFH and vertical 
profiles of each disc model (Czekaj et al. 2014| tab. 7 and 
fig. 4). The BesangonB model use the same local mass 
density as we do from JahreiB & Widen (1997) based on 
25 pc though they add the thick disc on top of this value 
which is an inconsistency as the JahreiB & Widen (1997) 
tabled accounts for all local stars not only the thin disc. 

Overall for the other disc mass models more mass seems 
to be sitting in stellar mass bins of bright stars which 
might be partly due to the local dwarf sample being almost 
2 times denser than the 200 and 100 pc sample at the 
upper MS (cf. table[^ Observations & CNS 5). Since we use 
volume complete samples to fit the luminosity function and 
take into account the scale height dilution according to the 
stellar ages we trust our mass distribution in the mass range 
that our observational evidence samples (0.5 — 10M©). 
Ideally we should have included a constraint for the low 
mass stellar mass density from other observations because 
now the mass factor and the low mass power law index 
are strongly correlated and certainly a bit too high. Better 
values for a more realistic two slope IMF would probably be 
in the lower left of the 1 a ellipse in figurej^ (e.g. mf ~ 1.05 
and ai cr 1.4). 


Table 7. Local mass density over stellar age in 10 "^Mq/pc^ 


Age [Gyr] 

Besangon A 

Besangon B 

JJ-model 

0 - 0.15 

20 

19 

9 

0.15 - 1 

55 

50 

36 

1 - 2 

46 

41 

30 

2 - 3 

33 

28 

27 

3 - 5 

58 

49 

52 

5 - 7 

61 

50 

54 

7 - 10 

117 

93 

84 

10 - 12 

- 

- 

46 

Thin disc 

390 

330 

338 

Thick disc 

29 

29 

17 

WD 

71 

71 

92t 

E 

490 

430 

447 


^ including brown dwarfs 


6 DISCUSSION 

We discuss our findings with respect to their model depen¬ 
dencies as well as our analysis method. Then we provide a 
comparison of the mass distribution from various IMFs and 
end this section with an empirically driven adaptation of our 
IMF to account for missing low mass star representation in 
our data. 


6.1 Mass factor 


The overall mass with our newly determined IMF compared 
to the JJ-model SFH increased by the thick disc fraction 
(6.5%) and the mass factor (mf = 1.09). The normalisation 
in paper I was done using ppdmf, total = 0.039 Mopc ^ from 
[JahreiB &: Widen (1997). Since we utilise new isochrones 
together with number densities derived from volume com¬ 
plete star counts a change of about 10 % is not unexpected 
but our value of 0.045 Mopc“^ is probably exaggerated. The 
present day mass fraction (getf) stays similar and also the 
remnant fraction is compatible with paper I values as shown 
in sect ion [53| though WD and faint stellar local mass density 
are arguably too high. The problem is that we do not have 
observational constraints for the whole mass range resulting 
in a degeneracy of the mass factor and the low mass slope 
(see figure]^ also valid for the high mass slope). We propose 
a solution to this in sectionlG^ 


6.2 Isochrones 


A crucial ingredient for our investigation is the used set of 
isochrones since it translates our analytical disc model into 
the realm of observables. When we were using the option 
provided originally by Galaxia ( Marigo et al.|2Q08 ) the high¬ 
est likelihood we could score was ln(£/Pmax) = —109 with 
slightly different IMF parameters. With the latest PARSEG 
isochrones the observations are much better fit by the model 
increasing the normed log-likelihood by 30. Still a few dis¬ 
crepancies are visible when inspecting figurej^ 

The cold dwarf V-Band problem was already mentioned in 
section 2.2 and is discussed in Ghen et aL] ( |2014| . We would 
not recommend to go much fainter in V band for such an 
analysis. 
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Section [5?T] discusses the huge likelihood penalty from the 
brightest magnitude bin in table[^ which could be due to 
unaccounted binaries in our sample, missing turn-off stars 
(which could be related to our discrete time steps), missing 
super solar metallicity stars or a high mass power law index 
which is too large. 

Tab le[^ shows missing stars in the 0th mag bin and an over 
abundance in the 1st for the giants. This is an indication for 
the red clump being too faint compared to the observations. 
Another striking feature is the complete absence of synthe¬ 
sised stars blueward from the red clump whereas Hipparcos 
shows several dozens. Maybe variable stars are not well rep¬ 
resented in our isochrones or it could again be a discrete 
time sampling effect. 

Of course not all differences are linked to the isochrones as 
our model SFH, AMR or IMF also affect the distribution of 
stars in the CMD. Another reason for mismatch is our used 
reddening law not accounting for inhomogeneous ISM which 
is probably the reason for all the unmatched faint giants in 
figure!^ 


6.3 Binarity 


In our observational sample we tried to account for all 
binaries that have Amag entries so that we could split them 
up. Apart from these another 365 stellar entries in our 
observational sample are listed in the Washington Double 


Star Catalogue (Mason et al. 2001) or the Catalogue of 


Components of Double and Multiple Stars (Dommanget & 
as multiple systems but are not resolved (i.e. 


they are not split up for this analysis but just kept as one 
star ). 

As mentioned in section [5T] we have a problem with binaries 
in massive stars which are hard to detect as lines are blurred 
and increased luminosity could also be due to ageing of a 
single star. Since binary fraction in massive stars is high 
and our likelihood gets a large penalty from the brightest 
dwarf mag bin not being matched well we believe that we 
miss quite a few high mass binaries. Other than that listings 


of binary stars in the Hipparcos catalogue (Perryman et al. 
1997) are said to be fairly eomplete ( [Lindegren et al.||I997 ) 


for AHp < 3.5 and an angular separation bigger than 
0.12 to 0.3arcsec (increasing with AHp). In CNS 5 all 
stellar systems from the literature with resolved magnitude 
differences are split into their respective components. 

The Besangon group is doing the favourable approach of 


accounting for binaries from the models side ( Czekaj et al. 
2014 Robin et al.|20I2] Arenou|20II ). They use an angular 


resolution limit of O.Sarcsec for resolved binaries in their 


Tycho 2 data. In Czekaj et al. (2014 fig. A.3) the relative 
difference in stellar mass frequency produced with binarity 
treatment in single, primary and secondary stars compared 
to the same sample excluding secondary stars is shown. 
With B-components included the IMF produces around 6 % 
more stars below I.IM© and around 6 % less above with a 
short transition in between. Overall the effect of binarity 
seems to play a secondary but not negligible role especially 
in massive dwarf stars. 


Table 8. Mass distribution, planetary nebula and supernova 
type II occurrences for different IMFs 


IMF: 

JJ 

Bes B 

KTG93 

Cha03 

Sa55 

Mass range in M0 


mass fraction 

in% 


8 - 100 

6 

4 

8 

22 

15 

1.4-8 

32 

30 

24 

29 

20 

I - 1.4 

13 

12 

9 

7 

5 

0.5 - I 

20 

23 

22 

16 

12 

0.08-0.5 

29 

31 

36 

26 

48 

PN* (1-8) 

994 

951 

707 

687 

467 

SNir (> 8) 

17 

12 

21 

47 

32 


^number of occurrences for an SSP with mass MiMF,25pc ~ 
4100 M0 at t = oo. 


6.4 Splitting the CMD 

Before exploiting the full 2d information of the CMD and 
adapting the statistical machinery as well as dealing with 
colour issues of the isochrones (not speaking of enhanced 
sensitivity to reddening) the easiest way to increase the data 
constraints is to split up the CMD into dwarf and giant 
sample. Weighting both into the final likelihood is a valuable 
gain since they represent stellar populations with different 
ages and masses (check the last columns of table[^&[^ and 
are still build from the same IMF (as well as SFH, AVR 
and AMR). Interestingly the likelihood penalty from both 
samples are similar though the dwarfs contribute 12 and 
the giants only 7 bins. This probably shows the challenge 
stars on the giant branch still pose to stellar evolutionary 
modelling. 

Exemplary for the insight gained from splitting up the CMD 
the 0 th mag bin can be inspected where a common sample 
would have balanced our model predicting too few giants 
and too many dwarfs. As these two bins have well distinct 
ages and masses constraints are put on different parts of the 
IMF. On the other hand the over abundance of stars in the 
3rd mag bin in giants and the 4th mag bin in both giants and 
dwarfs is indicating too many stars for the IMF around M©. 
This is in balance with the depletion of stars from around 
1.3 Mq (2nd mag in giants and 3rd in dwarfs) in order to 
fit the two slope power law with respect to our constructed 
likelihood. 


6.5 High mass slope 


Our high mass power law index of o = 3.02 is at the higher 
end of the literature values which can be partly attributed 
to the power law break being at a comparatively high mass. 
Still recent studies converge towards a steeper high mass 
slope as for example both new default IMFs of the Besangon 
model ( Czekaj et al.|2QT4 ) have similar values. 

An important measure for an IMF in terms of Galaxy sim¬ 
ulation is the fractional mass going into stars heavier than 
8 M 0 as they supposedly explode as supernova of typell 
(SNH) and have the highest and fastest stellar feedback to 
the Galactic evolution in terms of gross elemental synthesis 
and heating of the gas phase. In tablej^we list the fraction 
of total mass and total star count going into the SN H mass 
bin for all investigated IMFs and Salpeter with a = 2.3: The 
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fraction of mass going into stars ending their lives as SNII 
ranges from 4% for BesangonB ( Czekaj et al.|2QT4 ) to 22% 
for Chabrier (2003). 

In order to obtain sensible results modellers usually tune 
their feedback physics according to their used IMF. Chemi¬ 
cal enrichment models (OEMs) and their predicted element 
ratios are especially sensitive to the number of SNII events 
and the high mass slope. So using a Chabrier IMF will need 
substantially different stellar yields to obtain similar results 
(i.e. reproduce observations) compared to an CEM using the 
new BesangonB IMF. Same is true for cosmological hydro- 
dynamical simulations trying to reproduce the number of 
dwarf galaxies around spirals and so forth. 

Still our high mass slope is an extrapolation for stars heav¬ 
ier than about 10 M© as massive stars are too rare and not 
homogeneously distributed enough to be represented by our 
observational sample limited to 200 pc. The deficit of stars 
in our brightest dwarf bin might not only be due to missing 
binaries but also indicate a higher abundance of high mass 
stars than represented by our two slope IMF. 


6.6 Empirically motivated 3 slope IMF extension 

Our observational sample not constraining the high (volume 
too small) and low mass stars (sparse V Band data and no 
reliable isochrones) is deteriorating our inferred IMF param¬ 
eters. A good solution would be to include observationally 
based priors into the likelihood determination directing the 
MCMC simulation to a more physically motivated solution. 
Possible constraints directly connected to the shape of the 
IMF and the mass normalisation could be the SNII rate, 
HII regions, PN rate, WD number density, low mass stellar 
density. But it is not trivial to account for the uncertainty 
of those observations in the prior function. 

As an empirical fix we introduce a second power law break 
at 0.5 M© leaving the shape of the IMF above the same 
but changing it for lower masses. With a look at tablej^all 
bins should stay the same except for the 8 < My bin which 
should decrease to 0.017 Mq/pc^ being a new value derived 
from volume complete near infrared data of the CNS 5 which 
will soon be published (Just et al. 2015, in preparation). This 
deflates our mass model but still uses our high quality data 
for the other mass bins. The proposed IMF parameters ful¬ 
filling this additional constraint are: ao = 1.26, mo = 0.5, 
and mf = 1.03. Compared to our low mass slope of ai = 1.49 
the extension to lower masses is a bit shallower and quite 
similar to KTG 93 also having its low mass power law break 
at 0.5 M©. 


7 CONCLUSION 

We use solar neighbourhood stars to determine a new fidu¬ 
cial IMF within the framework of our local vertical Milky 
Way disc model (JJ-model). For that we carefully select vol¬ 
ume complete samples based on dereddened and binary cor¬ 
rected Hipparcos and CNS 5 data. Then we use Galaxia to 
create the corresponding mock observations from our JJ- 
model. We construct a likelihood by assuming a discrete 
Poisson process for the star count in magnitude bins dif¬ 
ferentiating between dwarfs and giants. With MCMC sim¬ 
ulations we sample the PDF of the two slope IMF param¬ 


eters. The derived IMF has a low mass power law index of 
ai = 1.49±0.08, a power law break at mi = 1.39±0.05 M©, 
a high mass index of a 2 = 3.02 ± 0.06, a mass factor of 
mf = 1.09 ± 0.04 with respect to our paper I mass nor¬ 
malisation. Except for physics not accurately represented 
in our model (binaries, inhomogeneous ISM, variable stars, 
low mass stellar atmospheres, metal enrichment law, thick 
disc) these findings are robust in our observationally backed 
mass range from 0.5 to 10 M©. An empirically driven low 
mass extension adds ao = 1.26 and mo = 0.5 and decreases 
the mass factor to 1.03. 

Independently from us the Besangon model using Tycho 2 
colour projections as observational constraints favours sim¬ 
ilar IMFs. The steep high mass slopes decrease the number 
of SNII ejected by an SSP compared to classical IMFs like 
Salpeter or Chabrier by a factor of about 3. 

The future of analytic Galaxy modelling will see increasing 
modularity to incorporate up-to-date theoretical progress in 
stellar atmospheres and evolutionary tracks. Observational 
biases like binarity, selection effects or reddening will be ac¬ 
counted for from the models side and the observational sam¬ 
ples will get diversified to overcome degeneracies in the var¬ 
ious model ingredients like SFH, AVR and CEM. With ever 
more realistic Milky Way models more elusive concepts like 
the dark halo, the spiral arm structure or chemical yields 
will get further constraints to fit within the holistic frame¬ 
work. 

To achieve that tools to measure the probability of model 
predictions given the data and schemes to optimise for the 
various data sets in an automated hierarchical fashion need 
to be implemented. 

In the future we plan to implement a detailed chemical en¬ 
richment scheme into our disc model and exploit the elemen¬ 
tal abundance information contained in surveys like RAVE, 
Gaia-ESO and SEGUE followed by the extension of the suc¬ 
cessful vertical description to different galactocentric radii. 
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